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The  concurrent  production  of  heat  and  electricity  within  residential  buildings  using  solid-oxide  fuel  cell 
(SOFC)  micro-cogeneration  devices  has  the  potential  to  reduce  primary  energy  consumption,  green¬ 
house  gas  emissions,  and  air  pollutants.  A  realistic  assessment  of  this  emerging  technology  requires 
the  accurate  simulation  of  the  thermal  and  electrical  production  of  SOFC  micro-cogeneration  devices 
concurrent  with  the  simulation  of  the  building,  its  occupants,  and  coupled  plant  components.  The  cali¬ 
bration  of  such  a  model  using  empirical  data  gathered  from  experiments  conducted  with  a  2.8  l<WAc  SOFC 
micro-cogeneration  device  is  demonstrated.  The  experimental  configuration,  types  of  instrumentation 
employed,  and  the  operating  scenarios  examined  are  treated.  The  propagation  of  measurement  uncer¬ 
tainty  into  the  derived  quantities  that  are  necessary  for  model  calibration  are  demonstrated  by  focusing 
upon  the  SOFC  micro-cogeneration  system’s  gas-to-water  heat  exchanger.  The  calibration  coefficients 
necessary  to  accurately  simulate  the  thermal  and  electrical  performance  of  this  prototype  device  are 
presented  and  the  types  of  analyses  enabled  to  study  the  potential  of  the  technology  are  demonstrated. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Micro-cogeneration  (also  known  as  residential  cogeneration 
and  small-scale  combined  heat  and  power)  is  the  concurrent  pro¬ 
duction  of  electricity  and  heat  from  a  single  fuel  source  with 
electrical  outputs  of  less  than  10-15  kW.  The  concurrent  production 
of  electrical  and  thermal  energy  has  the  potential  to  reduce  primary 
energy  consumption,  greenhouse  gas  emissions,  and  air  pollutants 
associated  with  providing  energy  services  to  buildings.  By  gener¬ 
ating  electrical  power  within  residential  buildings,  this  emerging 
technology  also  has  the  potential  to  reduce  electrical  transmission 
and  distribution  inefficiencies  and  to  alleviate  utility  peak  demand 
problems.  A  number  of  manufacturers  worldwide  are  developing 
micro-cogeneration  devices  based  upon  fuel  cells,  internal  com¬ 
bustion  engines,  and  Stirling  cycles  [1,2],  while  adoption  of  these 
technologies  is  being  encouraged  in  many  countries  through  finan¬ 
cial  incentives  and  favourable  electricity  tariff  structures. 

Micro-cogeneration  devices  have  only  modest  fuel-to-electrical 
conversion  efficiencies:  some  existing  prototypes  have  efficiencies 
as  low  as  5%  [3-5]  (net  AC  electrical  output  relative  to  the  source 
fuel’s  lower  heating  value  (LHV)).  Although  solid-oxide  fuel  cell 
(SOFC)  technologies  have  potential  to  deliver  electrical  efficien¬ 
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cies  as  high  as  45%  [6],  these  levels  have  not  yet  been  realized  in 
actual  prototype  devices.  These  electrical  efficiencies  are  relatively 
low  compared  to  combined-cycle  central  power  plants  (the  state- 
of-the-art  for  fossil-fuel-fired  central  power  generation),  which 
can  achieve  efficiencies  in  the  order  of  55%  [7,8].  Consequently,  it 
is  imperative  that  the  thermal  portion  of  the  micro-cogeneration 
device’s  output  be  well  utilized  to  supply  the  building’s  space  heat¬ 
ing,  space  cooling  (through  a  thermally  activated  cycle),  and/or 
domestic  hot  water  (DHW)  heating  needs.  If  this  thermal  out¬ 
put  cannot  be  exploited,  micro-cogeneration  will  not  compare 
favourably  to  the  best-available  central  power  generation  technolo¬ 
gies. 

However,  accurately  analyzing  the  utilization  of  a  micro¬ 
cogeneration  device’s  thermal  output  is  complicated  by  strong 
coupling  between  the  micro-cogeneration  unit,  other  thermal  plant 
components,  and  the  building’s  thermal  and  electrical  demands. 
This  complexity  can  be  illustrated  with  a  simple  example  that 
considers  a  micro-cogeneration  unit  that  is  configured  to  follow 
a  house’s  electrical  loads.  Lighting  and  appliance  demands  may 
peak  late  in  the  evening,  resulting  in  substantial  thermal  output 
from  the  micro-cogeneration  unit.  However,  there  may  be  little 
demand  for  space  heating  at  this  time  as  the  house  is  allowed  to  cool 
slightly  during  the  night.  Similarly  there  may  be  little  demand  for 
DHW.  Consequently,  the  system  will  likely  integrate  some  storage 
device  to  hold  the  thermal  energy  until  a  demand  exists.  The  vol¬ 
ume  and  thermal  characteristics  of  the  storage  tank,  the  occupant 
electrical  and  hot  water  usage  patterns,  the  house’s  thermal  charac- 
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teristics,  and  prevailing  weather  all  influence  whether  this  thermal 
energy  output  of  the  micro-cogeneration  device  will  be  exploited 
or  wasted.  The  potential  design  and  operational  combinations  of 
these  factors  are  almost  limitless.  In  order  to  accurately  assess  the 
performance  of  these  emerging  technologies,  therefore,  it  is  imper¬ 
ative  that  models  of  micro-cogeneration  devices  be  incorporated 
into  whole-building  simulation  tools  that  account  for  the  types  of 
interactions  outlined  above. 

These  factors  motivated  the  formation  of  Annex  42  of  the  Inter¬ 
national  Energy  Agency’s  Energy  Conservation  in  Buildings  and 
Community  Systems  Programme  (IEA/ECBCS)  [5].  The  objectives 
of  this  international  collaborative  effort  were  to  develop  simula¬ 
tion  models  that  advance  the  design,  operation,  and  analysis  of 
micro-cogeneration  systems,  and  to  apply  these  models  to  assess 
the  technical,  environmental,  and  economic  performance  of  the 
technologies.  These  objectives  were  accomplished  by  developing 
and  incorporating  models  of  micro-cogeneration  devices  within 
existing  whole-building  simulation  programs  such  as  ESP-r  [9], 
EnergyPlus  [10],  and  TRNSYS  [11  ].  These  models  are  more  detailed 
than  the  simple  performance-map  methods  that  have  been  pre¬ 
viously  applied  to  assess  micro-cogeneration  and  which  cannot 
accurately  treat  the  thermal  coupling  to  the  building  and  its  thermal 
plant  as  outlined  above  (e.g.  [12-24]). 

One  of  the  mathematical  models  developed  by  IEA/ECBCS  Annex 
42  is  suitable  for  simulating  the  performance  of  SOFC  micro¬ 
cogeneration  devices  [25,26].  This  is  a  system-level  model  that 
considers  the  thermodynamic  performance  of  all  components 
that  consume  energy  and  produce  the  SOFC  micro-cogeneration 
device’s  thermal  and  electrical  output.  The  model  requires  calibra¬ 
tion  using  empirical  data  that  can  be  acquired  from  the  testing  of 
coherent  systems  or  components  and  is  designed  for  operation  at 
a  time  resolution  that  is  in  the  order  of  minutes. 

This  paper  demonstrates  how  the  IEA/ECBCS  Annex  42  SOFC 
micro-cogeneration  model  can  be  calibrated  using  empirical  data. 
This  demonstration  is  performed  using  data  gathered  through 
experiments  conducted  on  a  prototype  2.8kWAC  SOFC  micro¬ 
cogeneration  system.  Pertinent  aspects  of  the  model  are  first 
reviewed.  The  experimental  equipment  and  measurement  meth¬ 
ods  employed  to  gather  the  calibration  data  are  then  presented. 
The  propagation  of  measurement  uncertainty  into  derived  quanti¬ 
ties  is  treated  and  the  methods  utilized  to  calibrate  model  inputs 
based  upon  these  data  are  then  elaborated.  The  calibration  param¬ 
eters  which  can  be  used  in  conjunction  with  the  model  to  simulate 
the  performance  of  this  specific  prototype  are  presented  and  then 
simulations  are  performed  to  demonstrate  the  use  of  the  calibrated 
model.  Concluding  remarks  are  then  provided. 

2.  Model  description 


The  system  will  include  a  power  conditioning  unit  to  convert  the 
SOFC’s  DC  electrical  output  to  AC.  All  SOFC  micro-cogeneration  sys¬ 
tems  will  include  a  heat  recovery  device  that  transfers  the  heat  of 
the  hot  product  gases  to  the  building’s  thermal  plant  to  provide 
space  heating,  space  cooling,  and/or  DHW  heating. 

As  previously  mentioned,  numerous  researchers  have  employed 
performance-map  methods  to  study  the  integration  of  SOFC  micro¬ 
cogeneration  systems  within  buildings  (e.g.  [12-15,20,22]).  These 
methods  employ  a  performance  map  (derived  either  from  empir¬ 
ical  evidence  or  from  detailed  modelling  performed  outside  the 
context  of  whole-building  simulation)  that  decouples  the  electrical 
and  thermal  performance  of  the  micro-cogeneration  device  from 
the  rest  of  the  thermodynamic  system.  In  contrast  to  this,  the  model 
that  is  the  subject  of  the  current  paper  explicitly  treats  the  inter¬ 
actions  between  the  micro-cogeneration  device  and  the  building 
and  its  occupants  on  a  time-step  basis  (typically  in  the  order  of 
minutes).  It  is  worth  noting  that  this  model  is  much  more  sim¬ 
plistic  than  many  detailed  SOFC  models  presented  in  the  literature 
(e.g.  [27,28]),  most  of  which  are  not  well  suited  for  the  purposes 
of  evaluating  the  energy  performance  of  SOFC  micro-cogeneration 
devices  since  they  focus  on  single  cells  or  stacks  of  cells  while  other 
components  (refer  to  Fig.  1 )  are  left  untreated. 

The  model  discretizes  the  SOFC  micro-cogeneration  system  into 
groupings  of  components  that  comprise  major  sub-systems,  such 
as  those  that  produce  electrical  power,  supply  air,  capture  heat  from 
hot  product  gases,  etc.  In  this  manner,  once  the  model  is  calibrated 
for  a  specific  SOFC  micro-cogeneration  device  analyses  can  be  con¬ 
ducted  to  explore  the  benefits  of  improving  the  performance  of 
individual  sub-systems.  For  example,  the  impact  of  improving  the 
heat  recovery  device  upon  overall  system  performance  can  be  simu¬ 
lated  without  recalibrating  the  portions  of  the  model  that  represent 
the  fuel  cell  power  module  (FCPM),  power  conditioner,  and  other 
sub-systems.  Additionally,  such  a  structure  facilitates  the  future 
development  of  more  detailed  modelling  methods  for  specific  sub¬ 
systems. 

The  model  discretizes  the  SOFC  micro-cogeneration  system  into 
nine  control  volumes.  Each  control  volume  is  modelled  in  as  rigor¬ 
ous  a  fashion  as  possible  given  the  constraints  of  computational 
efficiency  and  the  need  to  calibrate  model  inputs  based  upon  the 
testing  of  coherent  systems.  It  is  worth  noting  that  the  equations 
described  in  this  section  will  be  recalculated  hundreds  of  thousands 
of  times  to  perform  a  single  annual  simulation  at  2  min  time-steps. 

The  following  sub-sections  illustrate  the  modelling  methods 
employed  by  focusing  on  two  control  volumes.  One  represents  the 
exhaust-gas-to-water  heat  exchanger  that  produces  the  useful  ther¬ 
mal  output.  The  other  control  volume  represents  the  fuel  cell  power 
module  which  includes  the  stack,  after-burner,  and  the  other  com¬ 
ponents  enclosed  by  the  dashed  line  in  Fig.  1. 


This  section  describes  only  those  aspects  of  the  SOFC  micro¬ 
cogeneration  model  that  are  the  subject  of  the  calibration 
procedures  presented  in  Section  5.  The  interested  reader  is  referred 
to  [25]  and  [26]  for  a  comprehensive  exposition  of  the  model. 

2.1.  Model  topology 

Fig.  1  illustrates  one  possible  system  configuration  of  a  SOFC 
micro-cogeneration  device.  It  is  important  to  note  that  the  fuel  cell 
stack  is  only  a  single  component  within  a  complex  energy  conver¬ 
sion  system.  Besides  the  fuel  cell  stack  (shown  in  grey),  the  system 
might  include:  an  afterburner  to  combust  unreacted  fuel,  an  air  fil¬ 
ter  and  pre-heater,  and  a  fuel  desulfurizer,  pre-heater,  pre-reformer, 
and  reformer.  A  compressor  may  be  required  to  supply  pressurized 
fuel  while  a  blower  will  likely  be  present  to  supply  air  to  provide 
oxygen  to  support  the  electrochemical  and  combustion  reactions. 


2.2.  Exhaust-gas-to-water  heat  exchanger 


A  schematic  representation  of  the  control  volume  encapsulating 
the  device  that  transfers  heat  from  the  hot  gases  exiting  the  FCPM 
control  volume  to  the  water  loop  connected  to  the  building’s  plant 
is  shown  in  Fig.  2.  The  state  point  labels  shown  in  the  figure  are 
used  in  the  development  that  follows. 

The  sensible  component  of  the  heat  transfer  from  the  hot  gases 
to  the  water  is  characterized  with  the  log  mean  temperature  differ¬ 
ence  (LMTD)  method  for  counterflow  heat  exchangers, 


^sensible  —  (tJA)eff 


(7g-in  —  7w-out)  —  (T'g-out  —  T'w-in) 


In 


(Tg_[n— Tw-out  \ 
7g-out-7w-in  J 


(1) 


where  Tg_in  is  the  temperature  of  the  hot  gases  exiting  the  FCPM 
and  entering  the  heat  exchanger  and  Tg-0ut  is  the  temperature  of 
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Fig.  1.  One  possible  system  configuration  of  a  SOFC  micro-cogeneration  device. 


the  cooled  gases  exiting  the  heat  exchanger.  Tw_in  is  the  temper¬ 
ature  of  the  cold  water  at  the  heat  exchanger  inlet  and  Tw_out  is 
the  temperature  of  the  warmed  water  exiting  the  heat  exchanger. 
(LM)eff  is  the  effective  product  of  the  heat  transfer  coefficient  and 
area  (WIC1). 

If  it  is  assumed  that  heat  loss  from  the  heat  exchanger  to  the 
ambient  is  negligible  and  that  the  heat  capacity  of  each  fluid  stream 
remains  constant  through  the  heat  exchanger,  then  the  following 
energy  balance  can  be  written  for  the  heat  transfer  between  the 
fluid  streams, 

^sensible  —  (N^p)g-in(^g-in  —  "Eg- out)  =  (NCp)w-in(Ew-out  —  ^w-in) 

(2) 

where  cP  is  the  fluid’s  molar  heat  capacity  (Jkmol-1  K-1 )  and  N  is 
its  molar  flow  rate  (kmol  s-1 ). 


The  (UA)eff  value  required  by  Eq.  (1)  must  be  evaluated  at  each 
time-step  of  the  simulation.  An  approach  is  employed  which  casts 
it  as  a  parametric  relation  of  the  water  and  product  gas  flow  rates1 , 

(UA)eff  =  hxs  0  +  hxs  | Nw  +  hxsaNl,  +  hxs  3Ng  +  hxsANj  (3) 

where  hxsi  are  empirical  coefficients  characterizing  a  specific  heat 
exchanger.  The  form  of  Eq.  (3)  facilitates  the  determination  of  the 
hxSj  coefficients  from  experimental  data,  as  will  be  shown  in  Sec¬ 
tion  5. 

The  latent  component  of  the  heat  transfer  from  the  hot  gases  to 
the  water  is  given  by, 

^latent  —  ^H20-cond^H20,fg  (4) 


1  Alternate  approaches  are  available,  as  treated  in  [25]  and  [26]. 
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where  NH2o-cond  is  the  rate  of  condensation  of  water  from  the  gas 
stream  (kmol  s-1 )  and  fiH2o,fg  is  the  molar  heat  of  vapourization  of 
water  (J  kmol-1 ). 

The  rate  of  condensation  is  expressed  in  a  parametric  form  that 
facilitates  the  determination  of  its  coefficients  from  empirical  data. 
The  functional  form  of  this  parametric  equation  was  established 
by  recognizing  that,  for  a  given  heat  exchanger  design,  the  rate  of 
condensation  will  be  primarily  influenced  by  the  concentration  of 
water  vapour  in  the  gas  stream  and  by  the  difference  between  the 
heat  exchanger’s  temperature  and  the  gas’  dew  point, 


where  f/fue i  and  Hair  represent  the  enthalpy  carried  into  the  control 
volume  by  fuel  and  air  (for  electrochemical  and  combustion  reac¬ 
tions  as  well  as  excess  air).  Pel  is  the  net  DC  power  production,  that 
is  the  stack  power  less  ohmic  losses  in  cabling  and  the  power  draw 
of  ancillaries  (e.g.  the  fan  that  supplies  the  air).  Hg_ in  represents 
the  enthalpy  carried  out  of  the  control  volume  by  the  exiting  gas 
stream,  that  is  the  hot  products  of  the  electrochemical  and  combus¬ 
tion  reactions  as  well  as  the  excess  air  and  the  inert  constituents 
of  the  fuel.  (This  is  also  the  inlet  condition  for  the  heat  exchanger, 
as  indicated  in  Eq.  (1).)  The  final  two  terms  in  Eq.  (7)  represent 
thermal  losses:  qSkin-ioss  is  the  radiant  and  convective  heat  transfer 
to  the  containing  room  while  qFCPM-to-diiution  represents  the  heat 
transfer  from  the  FCPM  to  the  air  stream  which  is  drawn  through 
the  micro-cogeneration  device’s  cabinet  to  comply  with  gas  venting 
requirements  of  safety  codes. 

As  is  common  in  the  literature  (e.g.  [28,12])  the  FCPM’s  electrical 
efficiency  (eei)  is  expressed  as  the  DC  power  production  relative  to 
the  fuel’s  LHV, 


e  NfuelLHVfuei 

where  Nfuel  is  the  rate  of  fuel  supply  to  the  FCPM  (kmol  s-1 ). 

Since  the  model  groups  the  fuel  cell  stack  with  other  compo¬ 
nents  such  as  the  afterburner  and  heat  exchangers  into  the  FCPM 
control  volume,  it  makes  no  attempt  to  simulate  the  electrochem¬ 
ical  processes  occurring  within  the  fuel  cell,  but  rather  represents 
the  performance  of  the  FCPM  using  a  parametric  relation  between 
the  electrical  efficiency  and  the  net  electrical  power  output3, 

Se\  =  [£o+£\Pe\+S2Pl\\  (9) 

The  term  in  Eq.  (7)  for  the  enthalpy  of  the  inlet  air  stream  rep¬ 
resents  a  summation  of  the  enthalpies  of  the  constituent  gases, 


i^H20-cond  —  (^cond- threshold  ^w-in) 


hx  i 


nh20 


N, 


g-in 


+  fix, - 


Hh20 


N, 


g-m 


(5) 


where  hx\  ,•  are  empirical  coefficients  characterizing  a  specific  heat 
exchanger.  NH2o  in  Eq.  (5)  is  the  molar  flow  rate  of  water  vapour  in 
the  gas  stream  entering  the  heat  exchanger  and  Ng_in  is  the  molar 
flow  rate  of  all  constituents  of  the  gas. 

^cond-threshoid  is  a  user-specified  fixed  value  that  represents 
the  threshold  of  the  water-inlet  temperature  above  which  con¬ 
densation  will  not  occur.  When  Tw_ in  is  below  Tcond_threshold  the 
condensation  rate  will  be  determined  with  Eq.  (5).  And  when  Tw_in 
is  above  Tcond_threshold  it  is  assumed  that  no  condensation  occurs. 

The  total  heat  recovery  from  the  hot  gases  to  the  water  (qHx)  is 
determined  by  summing  the  results  of  Eqs.  (1)  and  (4). 


QUX  —  ^sensible  +  ^latent  (6) 

Section  5  details  the  methods  used  to  determine  the  hxs  l  coeffi¬ 
cients  for  Eq.  (3)  and  the  hxu  coefficients  and  Tcond_threshold  for  Eq. 
(5). 


2.3.  Fuel  cell  power  module 

The  energy  balance  for  the  FCPM  (indicated  by  the  dashed  line 
in  Fig.  1 )  can  be  written  in  the  following  form2, 

^fuel  +  ^air  —  ^el  +  f^g-in  +  Qsldn-loss  +  ^FCPM-to-dilution  (7) 


Hair  = 


(10) 


where  is  the  molar  enthalpy  (J  kmol-1 )  and  Nz  is  the  molar  flow 
rate  (kmol  s-1 )  of  gas  constituent  i  (e.g.  N2,  02). 

As  the  form  of  Eq.  (9)  includes  an  implicit  relationship  between 
the  fuel  and  air  supply  rates,  the  model  employs  an  explicit  rela¬ 
tionship  to  determine  the  FCPM’s  excess  air  ratio.  The  excess  air 
ratio  (A)  is  given  by, 


A.=  ^-l  (11) 

Ns. 

air 

where  Nair  is  the  air  supply  rate  of  the  FCPM  and  equals 
(kmols-1).  N^ir  is  the  air  supply  rate  required  to  support  stoi¬ 
chiometric  electrochemical  and  combustion  reactions  of  the  fuel 
constituents. 

A  parametric  function  is  used  to  relate  the  excess  air  ratio  to  the 
net  electrical  output4, 


^  _  (a0  +  a\Pe\  +  a2^eiX^  +  a3^air)  ^ 

~~  iV7 

air 

where  Tair  is  the  temperature  of  the  air  (°C)  supplied  to  the  FCPM. 

The  form  of  Eq.  (12)  facilitates  the  determination  of  the  a2  coef¬ 
ficients  from  experimental  data,  as  will  be  shown  in  Section  5. 

Eqs.  (10-12)  along  with  a  polynomial  equation  that  relates  the 
enthalpy  of  gas  constituents  to  temperature  are  used  to  estab¬ 
lish  the  Hair  term  of  Eq.  (7)  each  time-step  of  the  simulation.  The 


-  3  The  model  includes  terms  to  represent  the  operational  degradation  as  well  as 

2  Only  terms  that  are  relevant  to  the  particular  SOFC  device  treated  in  this  paper  the  impact  of  stop-start  cycling,  but  these  are  not  shown  here  for  the  sake  of  clarity, 
are  shown  here  for  the  sake  of  clarity.  4  Alternate  approaches  are  available,  as  treated  in  [25]  and  [26]. 
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Fig.  3.  Experimental  configuration  to  control  flow  rate  and  temperature  of  water  entering  the  heat  exchanger. 


remaining  terms  of  this  energy  balance  are  resolved  with  similar 
techniques,  which  is  then  re-arranged  and  solved  each  time-step  of 
the  simulation  to  predict  the  temperature  of  the  hot  gases  entering 
the  exhaust-gas-to-water  heat  exchanger  (refer  to  Tg_in  in  Eq.  (7)). 

The  next  section  describes  the  experiments  that  were  conducted 
with  a  prototype  2.8 1<WAC  SOFC  micro-cogeneration  device  in  order 
to  gather  the  data  required  to  calibrate  the  model’s  empirical  coef¬ 
ficients  such  as  those  of  Eqs.  (3),  (5),  (9),  and  (12). 


3.  Experimental  and  measurement  procedures 

An  experimental  programme  was  designed  to  examine  the 
performance  of  the  prototype  2.8  kWAC  SOFC  micro-cogeneration 
device  under  a  series  of  controlled  boundary  conditions  accord¬ 
ing  to  a  protocol  established  by  IEA/ECBCS  Annex  42  [29].  Many  of 
these  tests  required  control  over  the  water  flow  rate  through  the 
SOFC  micro-cogeneration  device’s  heat  exchanger  and  the  water 
inlet  temperature.  The  experimental  set-up  used  to  achieve  this 
control  is  shown  schematically  in  Fig.  3. 

Water  was  pumped  from  a  storage  tank  to  the  micro¬ 
cogeneration  device’s  heat  exchanger.  From  there  the  water  flowed 
through  a  fan-coil  before  returning  to  the  storage  tank.  As  the 
circulating  pump  was  operated  at  constant  speed,  the  flow  rate 
of  water  through  the  micro-cogeneration  device’s  heat  exchanger 


was  controlled  by  manually  setting  a  throttling  vale.  An  isola¬ 
tion  valve  downstream  of  the  pump  was  manually  controlled  to 
increase  back  pressure,  enabling  a  further  reduction  in  the  water 
flow  rate  through  the  heat  exchanger.  The  lowest  steady  water  flow 
rate  through  the  heat  exchanger  that  could  be  be  achieved  was 
4  L  min-1 .  The  highest  flow  rate  was  limited  by  the  pump’s  capacity 
and  was  approximately  equal  to  12  L  min-1. 

The  fan-coil  was  used  to  dissipate  heat  from  the  loop  when  the 
desired  water  temperature  was  greater  than  that  of  the  room  air. 
An  on-off  controller  with  a  0.2  °C  dead-band  cycled  the  fan-coil  on 
when  necessary  to  achieve  the  desired  water  inlet  temperature  at 
the  micro-cogeneration  device’s  heat  exchanger.  This  resulted  in 
a  small  degree  of  oscillation  although  the  control  was  mostly  sat¬ 
isfactory.  When  the  desired  water  temperature  was  below  that  of 
the  room  air,  warm  water  was  drained  downstream  of  the  micro¬ 
cogeneration  device.  This  volume  of  water  was  replenished  by 
adding  cold  water  from  the  mains  to  the  tank.  The  minimum  heat 
exchanger  water  inlet  temperature  was  thus  regulated  by  the  tem¬ 
perature  of  the  water  mains  (approximately  6°C).  The  maximum 
temperature  was  restricted  to  60  °C  in  order  to  protect  the  circulat¬ 
ing  pump. 

Once  steady  conditions  were  achieved,  measurements  were 
logged  to  file  for  a  period  of  time  to  provide  sufficient  data  to  ana¬ 
lyze  the  statistical  variation  of  the  measured  and  derived  quantities 
for  each  test.  Fig.  4  illustrates  the  flow  rate  and  heat  exchanger  water 
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Fig.  4.  Tw  _in  and  Nw  over  duration  of  one  test. 


inlet  temperature  for  the  test  that  was  configured  to  supply  30  °C 
water  to  the  heat  exchanger  at  the  lowest  flow  rate  possible.  As  can 
be  seen  from  the  graph,  ideal  steady  conditions  could  not  be  main¬ 
tained  over  the  duration  of  the  test.  Control  over  the  water  flow  rate 
was  found  to  be  more  stable  than  that  over  the  water  inlet  temper¬ 
ature.  In  general,  steady  thermal  conditions  were  more  difficult  to 
achieve  at  lower  entering  water  temperatures.  Notwithstanding, 
the  variations  in  the  water  inlet  temperature  were  deemed  to  be 
acceptable.  The  impact  of  these  variations  upon  derived  quantities 
will  be  examined  in  Section  4. 

In  addition  to  experiments  that  controlled  thermal  boundary 
conditions,  tests  were  performed  in  which  the  SOFC  micro¬ 
cogeneration  device’s  electrical  output  was  modulated.  This  was 
achieved  by  varying  the  stack  current  demanded  by  the  device’s 
internal  controller. 

The  micro-cogeneration  device  and  the  water  loop  were 
instrumented  to  record  both  electrical  and  thermal  conditions 
throughout  the  tests.  Voltage  and  current  were  measured  at  the 
points  where  power  flowed  to  the  power  conditioning  system,  to 
the  battery,  and  to  the  DC-powered  ancillary  devices  (refer  to  Fig.  1 ). 
The  AC  output  from  the  power  conditioning  system  was  also  instru¬ 
mented  as  were  the  AC-powered  ancillary  devices.  Voltage  taps 
were  placed  to  measure  DC  voltage  at  the  stack  exit  (i.e.  at  the  start 
of  the  transmission  cable  carrying  power  to  the  PCU)  and  at  the  AC 
ancillary  devices.  A  current  shunt  was  installed  to  measure  the  total 
ancillary  current  draw.  A  watt  transducer  was  used  to  monitor  the 
AC  output  to  the  grid. 

The  flow  rates  of  natural  gas  supplied  to  the  FCPM’s  stack  and 
burner  (fired  to  maintain  stack  temperatures  when  necessary)  were 
measured  independently  using  two  mass  flow  controllers.  Two  ven¬ 
turi  pressure  transducers  were  used  to  measure  the  flow  rates  of 
air  to  the  stack  and  burner. 

The  flow  rate  of  water  through  the  heat  exchanger  was  measured 
at  its  inlet  using  a  turbine  water  flow  meter.  Type-T  thermocouples 
were  used  to  measure  the  temperature  of  the  water  at  the  heat 
exchanger  inlet  and  outlet.  Gas  temperatures  were  measured  at  the 
heat  exchanger  inlet  and  outlet  using  type-K  thermocouples. 

Due  to  the  heat  exchanger’s  design,  when  water  vapour  con¬ 
densed  from  the  exhaust  gases  the  water  droplets  would  drip  onto 
the  thermocouple  measuring  Tg_in  (refer  to  Fig.  2).  This  resulted  in 


erroneous  temperature  readings,  a  fact  that  did  not  hinder  model 
calibration  efforts  but  rather  assisted  in  identifying  Tcond_threshold 
in  Eq.  (5),  as  will  be  treated  in  Section  5.  The  micro-cogeneration 
device  collects  the  condensate  in  an  internal  reservoir.  When  full,  a 
float  valve  triggers  a  pump  to  drain  this  reservoir.  A  rain  gauge  tilt 
bucket  was  located  to  collect  the  pumped  condensate  to  measure 
its  volumetric  flow  rate.  This  gauge  was  calibrated  to  tilt  for  each 
accumulation  of  8.24  mL  of  condensate. 

The  micro-cogeneration  device  is  designed  such  that  the  cooled 
gases  exiting  the  heat  exchanger  are  mixed  with  the  dilution  air 
that  is  drawn  through  the  cabinet  to  control  skin  losses  to  the 
containing  room.  The  temperature,  velocity,  and  relative  humid¬ 
ity  of  these  mixed  gases  were  measured  downstream  of  the  mixing 
point.  A  velocity  probe  was  used  to  measure  the  velocity  of  this  gas 
stream.  Due  to  the  configuration  of  the  micro-cogeneration  device’s 
exhaust  chimney  it  was  not  possible  to  take  these  measurements 
in  a  region  of  fully  developed  flow.  Rather,  measurements  had  to  be 
taken  close  to  a  90°  bend  in  the  duct.  During  the  exploratory  phase 
of  the  work,  the  probe  was  inserted  at  numerous  locations  across 
the  duct  and  the  measured  velocity  profile  examined  to  choose  the 
most  representative  location  to  mount  the  probe.  These  limitations 
resulted  in  significant  uncertainty  in  the  measured  flow  rate  of  the 
combined  gas  stream. 

Finally,  the  ambient  temperature  and  relative  humidity  in  the 
test  room  were  measured  approximately  1  m  above  the  top  of  the 
fuel  cell  enclosure  and  approximately  1  m  away  from  the  air  inlet 
side  of  the  micro-cogeneration  device. 

As  detailed  in  the  next  section,  the  bias  and  precision  errors  from 
the  primary  measurements  outlined  above  (e.g.  temperatures,  flow 
rates)  propagate  through  into  the  derived  quantities  (e.g.  ( LM)eff  of 
Eq.  ( 1 )).  In  order  to  minimize  the  bias  errors,  a  number  of  the  instru¬ 
ments  described  above  were  calibrated.  These  include  the  water 
flow  meter,  the  thermocouples  at  the  heat  exchanger’s  water  inlet 
and  outlet,  the  AC  power  flow  meter,  and  the  natural  gas  flow  meter. 
These  calibrations  were  effected  by  comparing  instrument  read¬ 
ings  to  reference  instruments  and  then  adjusting  offset  and  slope 
parameters  to  adjust  the  translation  of  voltage  signals  to  measured 
quantities. 

Instantaneous  measurements  of  the  FCPM’s  DC  power  produc¬ 
tion,  the  FCPM’s  air  and  fuel  supply  rates,  and  the  power  flow  to  the 
battery  were  taken  every  second  and  the  averages  over  the  minute 
were  logged  to  file.  All  other  measurements  were  taken  every  15  s 
and  the  four  values  averaged  to  log  the  data  at  each  minute.  The 
condensate  flow  rate  was  logged  at  the  same  frequency,  but  using 
a  separate  data  acquisition  system.  Each  of  these  measurements 
records  the  number  of  times  the  bucket  had  been  tilted  during 
the  preceding  minute.  The  time  stamps  in  each  file  were  used  to 
synchronize  the  measurements. 

Infrared  images  of  the  micro-cogeneration  device  were  captured 
during  one  test  at  which  the  device  was  producing  its  maximum 
power.  Three  of  the  four  side  faces  and  the  top  of  the  SOFC  micro¬ 
cogeneration  device’s  enclosure  provided  unobstructed  views  for 
the  imaging.  These  images  were  used  to  derive  thermal  contour 
maps  by  taking  into  account  the  surface  emissivities. 

A  gas  chromatograph  was  used  to  analyze  the  content  of  the 
natural  gas  supply  a  few  days  prior  to  the  experiments.  This  deter¬ 
mined  the  molar  fractions  of  each  constituent  of  the  gas  supply  in 
order  to  accurately  determine  its  lower  heating  value. 


4.  Calculation  of  calibration  data  and  analysis  of 
uncertainties 

The  previous  section  detailed  the  measurements  taken  during 
the  experiments.  In  order  to  calibrate  the  model,  these  primary 
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measurements  were  used  to  derive  the  variables  of  interest  to  the 
model.  This  section  details  the  calculation  of  these  derived  quanti¬ 
ties  and  their  associated  uncertainties.  The  methods  illustrated  here 
for  treating  the  exhaust-gas-to-water  heat  exchanger  equally  apply 
to  other  aspects  of  the  model,  such  as  the  electrical  efficiency  of  the 
FCPM,  the  DC-AC  conversion  efficiency  of  the  power  conditioning 
system,  etc. 

As  discussed  in  Section  2.2,  the  equation  relating  the  effective 
product  of  the  heat  transfer  coefficient  and  area  to  the  flow  rates  of 
water  and  gas  through  the  heat  exchanger  (refer  to  Eq.  (3))  must  be 
calibrated  from  the  experimental  data.  Referring  to  Eqs.  (1 )  and  (2), 
it  can  be  shown  that  (LM)eff  can  be  derived  from  five  of  the  primary 
measurements  described  in  Section  3, 


(LM)eff  = 


(Ncp)w_jn(7w 

-out  ^w-in) 


(Tg_in — Tw-out )— (rg_out — Tw_in ) 
ln(^  Vin^w-out\ 

y  'g-out-'w-in  J 


(13) 


Eq.  (13)  was  evaluated  for  each  minute  of  recorded  data  using 
the  four  temperature  readings  (Tw_ in,  Tw_ out,  Tg_in,  Tg_0 ut)  and  the 
water  flow  rate  measurement  (Nw_in).  The  heat  capacity  of  the 
water  entering  the  heat  exchanger  (cPiW-in)  was  calculated  from 
Tw_in  using  a  parametric  relation  [26]. 

The  method  recommended  by  the  American  Society  of  Mechani¬ 
cal  Engineers  (described  in  [30]  and  [31  ])  was  used  to  calculate  the 
uncertainties  of  the  measured  quantities  and  to  propagate  these 
uncertainties  into  the  derived  quantities.  With  this  a  bias  error 
was  assigned  to  each  primary  measurement.  These  were  estab¬ 
lished  based  upon  the  instrumentation  specifications  as  either  an 
absolute  error  as  a  percent  of  full-scale  measurement  and/or  a 
reading  error  as  a  percent  of  the  value  measured.  Where  instru¬ 
ments  were  calibrated  (refer  to  Section  3)  the  bias  error  was 
established  based  upon  the  calibration  parameters.  In  these  cases, 
the  bias  error  was  set  based  upon  either  the  average  or  maxi¬ 
mum  deviation  of  the  corrected  measured  values  to  the  reference 
values. 

In  some  cases  additional  bias  errors  were  assigned  based  upon 
judgement.  For  example,  a  substantial  bias  error  was  assigned  to 
the  velocity  measurement  of  the  combined  exhaust  gas  stream  due 
to  the  restrictions  on  instrument  placement,  as  discussed  in  Section 
3.  As  another  example,  an  additional  bias  error  was  assigned  to  the 
condensate  flow  rate  measurement.  As  described  in  Section  3,  con¬ 
densate  is  measured  by  a  rain  gauge  tilt  bucket  after  it  is  pumped 
from  an  internal  reservoir.  The  time  lag  between  the  pumping  and 
measurement  actions  introduced  some  uncertainty  to  the  conden¬ 
sate  flow  rate  measurement.  Consequently  a  bias  error  of  50  mL 
(the  approximate  volume  of  the  internal  reservoir)  was  assigned  to 
the  measurement  of  the  condensate  flow  over  the  duration  of  each 
experiment. 

The  total  bias  for  each  measurement  point  is  calculated  from 
the  individual  bias  errors  for  that  sensor  using  the  root-sum-square 
method, 


of  logged  readings, 


5  = 


Ef=i^-*)2 

N  —  1 


1/2 


(15) 


where  N  is  the  number  of  logged  readings.  It  is  worth  noting  that 
the  data  were  logged  at  1  -min  intervals  based  upon  either  1  s  or  15  s 
instantaneous  readings  (refer  to  Section  3).  TheX,  values  ofEq.  (15) 
are  the  1-min  averaged  values  since  the  instantaneous  data  were 
not  logged.  It  is  also  worth  noting  that  S  has  the  same  value  for  each 
data  point  within  a  given  test. 

The  precision  index  of  the  average  value  of  a  parameter  for 
a  given  test  is  lower  than  that  for  the  individual  measurements 
according  to, 


(16) 


Finally,  the  bias  and  precision  indices  are  combined  to  express 
the  uncertainty  in  a  measured  quantity  at  the  95%  confidence  level, 


U95%  =  V/B2+(  tS)2 


(17) 


where  t  is  the  standard  statistical  Student’s  t-value  and  is  a  function 
of  the  value  of  N  used  in  evaluating  Eq.  (15). 

The  uncertainty  of  a  derived  quantity  is  determined  by  propa¬ 
gating  the  bias  and  precision  indices  of  the  measurements  that  are 
used  to  calculate  the  derived  quantity.  For  example,  the  bias  error 
for  (UA)eff  is  calculated  as  follows  (refer  to  Eq.  (13)), 
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(18) 


The  precision  index  for  (UA)eff  is  determined  in  a  similar  manner 
and  the  overall  uncertainty  determined  using  Eq.  (17). 

The  propagation  of  measurement  uncertainties  into  Eq.  (18)  is 
demonstrated  by  examining  the  test  that  was  illustrated  in  Fig.  4. 
The  bias  errors  and  precision  indices  for  the  four  temperature 
and  one  water  flow  rate  measurements  used  in  the  equation  are 
summarized  in  Table  1.  The  bias  errors  reported  in  the  table  are 
the  average  for  the  82  measurement  points  of  the  test.  Likewise, 
the  precision  index  is  that  corresponding  to  each  individual  mea¬ 
surement,  and  not  the  precision  index  of  the  set  average  (i.e.  it 
represents  S  of  Eq.  (15),  not  5avg  ofEq.  (16)). 

Table  1  also  lists  the  average  uncertainty  at  the  95%  confidence 
level  of  the  82  measurements  of  each  of  the  five  parameters.  As  can 
be  seen,  the  precision  indices  are  the  predominant  determinant  of 
the  uncertainty  of  the  water  temperature  measurements,  an  obser- 


B  =  [B21+Bl+--  +  B2k]V2  (14) 

For  each  of  the  tests  examining  the  performance  of  the  gas-to- 
water  heat  exchanger  the  desired  boundary  conditions  (e.g.  Tw_in 
and  Nw-in)  were  held  as  constant  as  possible  for  a  period  of  time 
and  data  logged  each  minute.  The  precision  index  of  a  single  mea¬ 
surement  within  a  given  test  is  calculated  based  on  the  average 
value  of  the  observed  parameter  during  that  test  and  the  number 


Table  1 

Uncertainty  parameters  for  test  at  Tw_in  =  30  °C  and  Nw  =  0.004  kmol/s. 


Measurement 

Average  value 
over  test 

B 

S 

^95% 

T/v-in 

30.60  °C 

0.10°C 

0.58  °C 

1.17  °C 

T/V-OUt 

43.38  °C 

0.10°C 

0.48  °C 

0.97  °C 

T  ■ 

1  g-m 

284.27  °C 

2.20  °C 

0.57  °C 

2.48  °C 

Tg-out 

45.04  °C 

2.20  °C 

0.36°C 

2.32  °C 

Nw  (kmols) 

4.0  x  10~3 

7.9  x  10“5 

2.7  x  10~5 

9.3  x  10“5 
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Fig.  5.  Derived  (LM)eff  values  and  associated  95%  error  bars  for  Tw_in  =  30 °C  and 
Nw  =  0.004  kmol/s. 

vation  consistent  with  the  Tw_in  measurements  plotted  in  Fig.  4.  In 
contrast,  the  instrument  bias  errors  are  the  predominant  determi¬ 
nants  of  the  uncertainty  of  the  gas  temperatures  and  the  water  flow 
rate. 

Eq.  (13)  was  applied  to  calculate  the  (LM)e ff  value  for  each  of  the 
82  measurement  points  of  the  test.  The  procedure  outlined  in  Eqs. 
(14),  (15),  (17),  and  (18)  was  then  applied  to  calculate  the  uncer¬ 
tainty  for  each  of  these  82  derived  (LM)e ff  values.  Fig.  5  plots  these 
derived  values  and  their  uncertainties.  The  test-averaged  (LM)e ff 
value  determined  from  the  82  measurement  points  and  its  error  bar 
are  also  shown  in  the  figure.  The  uncertainty  of  the  test-averaged 
(LM)eff  value  is  less  than  that  for  individual  measurements  due  to 

Eq*  (16). 

The  procedure  outlined  in  this  section  was  applied  to  each  test 
to  produce  sets  of  test-averaged  derived  quantities  (e.g.  (HA)e ff,  £el, 
X)  for  various  combinations  of  the  boundary  conditions  (e.g.  Tw_in, 
Nw,  Pel).  The  calibration  of  the  model  using  these  sets  of  data  is  the 
subject  of  the  next  section. 

5.  Model  calibration 

This  section  treats  the  calibration  of  the  model  inputs  described 
in  Section  2  using  the  measured  and  derived  data  described  in 
Sections  3  and  4. 


5.1.  Exhaust-gas-to-water  heat  exchanger 

The  experiments  yielded  17  derived  values  of  (HA)eff  at  vari¬ 
ous  water  (Nw)  and  gas  (Ng)  flow  rates.  A  non-linear  regression 
method  was  used  to  establish  the  values  of  the  hxs  i  coefficients 
that  produced  the  best  fit  of  these  data  to  Eq.  (3).  The  values  of  the 
coefficients  determined  from  this  analysis  are  presented  in  Table  2. 

It  is  important  to  note  that  the  test-averaged  derived  (HA )eff,  Nw, 
and  Ng  values  (refer  to  Section  4)  were  used  to  regress  these  coeffi¬ 
cients.  As  a  non-linear  regression  method  was  used  it  is  not  possible 
to  propagate  uncertainties  in  these  values  into  the  hxsi  coefficients 
using  the  methods  elaborated  in  Section  4.  However,  despite  the 
deterministic  appearance  of  the  coefficients  in  Table  2,  the  regres¬ 
sions  inherently  contain  uncertainty  and  this  level  of  uncertainty 
is  commensurate  with  the  uncertainty  of  the  derived  (HA)eff,  Nw, 
and  Ng  values. 

Fig.  6  compares  the  (HA)eff  determined  with  Eq.  (3)  and  the 
coefficients  of  Table  2  with  the  (HA)eff  values  derived  from  the 
measurements.  The  uncertainty  bars  determined  in  Section  4  are 
plotted  in  the  figure.  The  left  side  of  the  figure  provides  a  view  nor¬ 
mal  to  the  Ng  axis  while  the  right  side  provides  a  view  normal  to  the 
Nw  axis.  As  can  be  seen,  the  functional  form  of  Eq.  (3)  well  repre¬ 
sents  the  dependency  of  ( HA)eff  on  the  two  flow  rates.  The  calibrated 
values  lie  within  the  error  bars  for  each  of  the  17  data  points. 

Fig.  7  provides  another  indication  of  the  goodness  of  fit  between 
the  calibrated  (HA )eff  values  and  those  derived  from  measurements. 
The  coefficient  of  determination  (r2  value)  was  0.98.  The  aver¬ 
age  error  (difference  between  the  calibrated  (HA)eff  value  and  that 
derived  from  measurements)  was  1.9%  while  the  root-mean-square 
error  was  2.1%.  The  maximum  error  for  a  single  point  was  3.2%. 

A  number  of  tests,  in  addition  to  the  17  described  above,  were 
conducted  to  explore  the  operation  of  the  heat  exchanger  under 
condensing  conditions.  One  of  these  tests  was  configured  to  identify 
Econd-threshoid  of  Eq.  (5).  This  variable  represents  the  threshold  of  the 
water  inlet  temperature  above  which  condensation  does  not  occur. 
The  examination  of  the  tilt  bucket  readings  during  preliminary  test¬ 
ing  indicated  an  approximate  range  within  which  Tcond_threshold  lay. 
However,  each  of  these  tests  was  time  consuming.  As  elaborated 
in  Section  3,  the  tilt  bucket  instrument  was  filled  only  after  the 
cogeneration  device’s  internal  condensate  reservoir  became  filled 
and  was  pumped  out.  Steady  conditions  had  to  be  held  for  long  peri¬ 
ods  of  time  (in  some  cases  many  hours)  in  order  to  register  readings 
at  the  tilt  bucket. 

Section  3  explained  how  the  formation  of  condensation  from 
the  exhaust  gases  led  to  erroneous  Tg_in  thermocouple  readings. 
Advantage  was  taken  of  this  fact  to  calibrate  Tcond-thresh0id  f°r  Eq. 
(5).  By  controlling  the  water  loop  illustrated  in  Fig.  3,  Tw_in  was 
slowly  reduced  from  33  °C,  which  the  preliminary  testing  had  indi¬ 
cated  was  above  Tcond-threshoid.  For  the  FCPM’s  electrical  output 
exercised  in  this  test,  Tg_in  was  approximately  280  °C.  As  shown  in 
Fig.  8  the  thermocouple  produced  reliable  data  for  the  first  portion 


Table  2 

Calibrated  coefficients. 

Gas-to-water  heat  exchanger  sensible  coefficients  for  Eq.  (3) 

Gas-to- water  heat  exchanger  latent  coefficients  for  Eq.  (5) 

FCPM  electrical  efficiency  coefficients  for  Eq.  (9) 

FCPM  air  supply  coefficients  for  Eq.  (12) 

hxSt o  =  83.1 ;  hxsi  =  4  798;  hxs>2  =  -138  x  103 ;  hxSt3  =  -353.8  x  103 ;  hxsA  =  -5.15  x  108 
hx it  =  —1.96  x  10  4 ;  hx\2  =3.1  x  10  3;  Tconci_threshold  =  35 °C 
e0  =  0.642388;  e3  =  -1.619  x  10~4;  e2  =2.26007  x  10~8 

a0  =  1.50976  x  10~3 ;  a ^  =  -7.76656  x  lO”7 ;  a2  =  1.30317  x  10“10 ;  a3  =  2.83507  x  10~3 

Range  of  applicability: 

5.2  °C  <  Twater  jn  <  59.8  °C 

3.63  x  1CT3  kmol/s  <  Nwater  <  1.09  x  10~2  kmol/s 

3010  W<  Pd  <  3728  W 

10.0  °C<Tair<  24.2  °C 

Froom  ~  22  °C 
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of  the  test.  (It  should  have  read  approximately  280  °C  throughout 
the  test.)  The  thermocouple,  however,  began  producing  unreliable 
readings  once  Tw_in  was  reduced  to  23  °C.  This  event  indicated  the 
first  formation  of  liquid  water  which  dripped  onto  the  thermocou¬ 
ple.  Even  as  the  inlet  water  temperature  was  warmed  to  35  °C,  the 
thermocouple  continued  to  produce  unreliable  readings,  indicating 
that  condensation  continued  to  occur.  It  took  considerable  time  for 
the  thermocouple  readings  to  stabilize.  This  tends  to  indicate  that 
although  the  onset  of  condensation  requires  a  low  value  of  Tw_ in, 
once  condensing  conditions  have  been  achieved  condensation  can 
occur  at  warmer  temperatures.  Based  upon  this  test  it  was  decided 
to  set  Tconcj_tiiresi1oid  to  35  °C. 


A  series  of  seven  tests  were  then  conducted  at  various  water  flow 
rates  and  values  of  Tw_ in  in  order  to  establish  the  hx\  ,■  coefficients 
of  Eq.  (5).  Sufficient  time  was  allowed  in  each  test  to  achieve  steady 
conditions.  Due  to  practical  constraints,  however,  these  tests  could 
only  be  conducted  with  a  nearly  constant  ratio  of  water  vapour 
in  the  exhaust  gas  stream  (refer  to  NH2o/Ng-in  in  Eq.  (5))-  (This  is 
determined  by  the  FCPM’s  operating  point.) 

A  non-linear  regression  method  was  used  to  establish  the  val¬ 
ues  of  the  hx\  i  coefficients  that  produced  the  best  fit  to  Eq.  (5). 
As  elaborated  above,  Tcond_threshoid  was  set  to  35  °C  to  perform  this 
regression.  The  values  of  the  coefficients  determined  from  this  anal¬ 
ysis  are  presented  in  Table  2. 

Fig.  9  compares  the  NH2o-cond  determined  with  Eq.  (5)  and 
the  coefficients  of  Table  2  with  the  Nn2o-cond  values  derived  from 


time  (hours) 


Fig.  7.  Goodness  of  fit  between  calibrated  and  measured  (UA)eff. 


Fig.  8.  Identification  of  rcond_threshold. 
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Fig.  10.  Calibrated  £ei  versus  measurements. 


the  measurements.  The  coefficient  of  determination  (r2  value) 
was  0.96.  The  average  error  (difference  between  calibrated  val¬ 
ues  and  those  derived  from  measurements)  was  10-6kmols-15 
while  the  root-mean-square  error  was  1.2  x  10-6  kmols-1.  The 
maximum  error  for  a  single  point  was  2.1  x  10-6  kmols-1.  The 
uncertainty  bars  determined  in  Section  4  are  plotted  in  the  figure. 
As  can  be  seen,  the  functional  form  of  Eq.  (5)  reasonably  repre¬ 
sents  the  dependency  of  N^o-cond  upon  Tw_in.  The  calibrated  values 
lie  within  the  error  bars  for  five  of  the  seven  data  points.  The 
greatest  deviation  between  measurement  and  calibration  occurs 
at  7w-in  ^  30  °C  where  the  condensation  flow  rate  is  very  small. 

5.2.  Fuel  cell  power  module  electrical  efficiency 

Data  from  seven  experiments  were  used  to  calibrate  the  FCPM 
electrical  efficiency  to  the  net  DC  power  production  (refer  to  Eq. 
(9)).  A  non-linear  regression  method  was  used  to  establish  the  val¬ 
ues  of  the  Si  coefficients  that  produced  the  best  fit  to  Eq.  (9).  The 
values  of  the  coefficients  determined  from  this  analysis  are  pre¬ 
sented  in  Table  2.  This  table  also  presents  the  range  of  Pe i  for  which 
the  correlation  is  applicable. 

Fig.  10  compares  the  £ei  determined  with  Eq.  (8)  and  the  coeffi¬ 
cients  of  Table  2  with  the  £el  values  derived  from  the  measurements. 
The  D95%  uncertainty  bars  are  plotted  in  the  figure.  The  coefficient 
of  determination  (r2  value)  was  0.58.  The  average  error  (difference 
between  the  calibrated  £el  value  and  that  derived  from  measure¬ 
ments)  was  0.4%  (in  relative  terms)  while  the  root-mean-square 
error  was  0.6%.  The  maximum  error  for  a  single  point  was  1.2%. 

5.3.  Fuel  cell  power  module  air  supply 

Data  from  28  experiments  were  used  to  calibrate  the  excess  air 
ratio  of  the  FCPM  air  supply.  A  non-linear  regression  method  was 
used  to  establish  the  values  of  the  az  coefficients  that  produced  the 


5  To  place  these  numbers  in  context,  a  condensation  rate  of  10-6  kmol/s  results  in 
approximately  40  W  of  heat  transfer. 


best  fit  to  Eq.  (12).  The  values  of  the  coefficients  determined  from 
this  analysis  are  presented  in  Table  2.  This  table  also  presents  the 
range  of  Pel  and  Tair  for  which  the  correlation  is  applicable. 

Fig.  11  compares  the  Nair  determined  with  Eq.  (12)  and  the 
coefficients  of  Table  2  with  the  Nair  values  derived  from  the  mea¬ 
surements.  The  t/95%  uncertainty  bars  are  plotted  in  the  figure. 
The  coefficient  of  determination  (r2  value)  was  0.77.  The  aver¬ 
age  error  (difference  between  the  calibrated  Nair  value  and  that 
derived  from  measurements)  was  2.3%  while  the  root-mean-square 
error  was  2.8%.  The  maximum  error  for  a  single  point  was  5.6%.  As 
Fig.  11  illustrates,  for  some  of  the  28  points  the  deviation  between 


Fig.  11.  Calibrated  Nair  versus  measurements. 
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SOFC  micro-cogeneration  device 


the  calibration  and  the  experimentally  derived  values  exceeds  the 
measurement  uncertainty.  This  result  indicates  the  potential  for 
improvements  to  the  form  of  Eq.  (12),  although  additional  data 
under  a  greater  range  of  operating  conditions  would  be  required 
to  explore  this. 

6.  Demonstrated  application  of  calibrated  model 

The  process  elaborated  in  the  previous  section  was  repeated 
for  all  aspects  of  the  model.  The  interested  reader  is  referred  to 
[29]  which  provides  the  complete  set  of  calibration  parameters  to 
accurately  represent  the  performance  of  this  2.8 1<WAC  SOFC  micro¬ 
cogeneration  device. 

The  calibrated  model  can  now  be  used  to  accurately  simulate  the 
thermal  and  electrical  performance  of  this  prototype  device  when 
coupled  to  a  house’s  thermal  plant.  This  is  demonstrated  through 
the  use  of  the  ESP-r  simulation  program  [9],  which  includes  the 
IEA/ECBCS  Annex  42  model,  in  conjunction  with  the  calibration 
parameters  presented  in  this  paper.  An  ESP-r  simulation  model 
was  configured  to  represent  a  typical  20-30-  year-old  detached 
two-storey  Canadian  house  with  160  m2  of  conditioned  floor  area. 
When  conditioned  to  21  °C  and  simulated  with  Ottawa  weather 
data,  the  house  has  an  annual  space  heating  demand  of  78  GJ, 
which  is  typical  of  this  vintage  of  construction,  house  size,  and 
location. 

The  SOFC  micro-cogeneration  device  is  coupled  to  supporting 
plant  components  as  illustrated  in  Fig.  12  in  order  to  provide  the 
house’s  space  heating.  The  thermal  output  of  the  SOFC  micro¬ 
cogeneration  device  is  transferred  to  a  water  storage  tank  by  a 
circulating  pump.  When  the  house’s  temperature  drops  below 
20.5  °C,  the  second  pump  cycles  on  to  circulate  hot  water  from  the 
tank  through  a  water-to-air  heat  exchanger.  Likewise,  the  fan  cycles 
on  to  circulate  air  from  the  house  through  the  heat  exchanger  and 
delivers  this  warmed  air  to  heat  the  house. 

The  tank’s  make-up  burner  cycles  on  when  then  tank  temper¬ 
ature  drops  below  50  °C.  To  avoid  excessive  temperatures,  energy 


is  extracted  from  the  tank  and  rejected  to  the  environment  when 
its  temperature  rises  to  92  °C.  The  SOFC  micro-cogeneration  device 
is  assumed  to  operate  continuously  throughout  the  year  producing 
2.8  kW  of  AC  power.  It  is  assumed  that  the  house  can  export  power 
to  the  central  electrical  grid  when  there  is  a  surplus  of  power,  and 
import  from  the  grid  when  there  is  a  deficit. 

An  annual  simulation  conducted  with  this  configuration  using 
Ottawa  weather  data  predicted  that  9366  m3  of  natural  gas 
would  be  consumed  by  the  SOFC  while  another  745  m3  would  be 
consumed  by  the  make-up  burner.  Over  the  year,  the  SOFC  micro¬ 
cogeneration  device  would  generate  more  than  two  and  a  half 
times  the  electrical  demands  of  the  occupants  and  plant  equipment 
(pumps  and  fan),  resulting  in  a  net  export  of  53.4  GJ  of  electricity 
to  the  central  grid.  Due  to  temporal  mismatches  between  electri¬ 
cal  production  and  demand,  some  importation  from  the  grid  would 
also  be  required.  The  monthly  integrated  electrical  balance  is  illus¬ 
trated  in  Fig.  13. 

Some  important  observations  regarding  the  suitability  of  this 
configuration  can  be  drawn  by  examining  the  monthly  integrated 
thermal  energy  balance  on  the  water  tank  (see  Fig.  14).  As  can 
be  seen,  during  the  winter  months  the  SOFC  micro-cogeneration 
device’s  thermal  output  is  insufficient  to  meet  the  house’s  space 
heating  demand.  Consequently,  extensive  use  of  the  tank’s  make-up 
burner  is  required.  During  this  period,  the  tank  never  overheats  and 
thus  no  dumping  of  heat  is  required.  During  the  summer,  however, 
the  vast  majority  of  the  SOFC  micro-cogeneration  device’s  thermal 
output  is  wasted.  There  is  also  significant  dumping  of  heat  during 
the  spring  and  autumn  due  to  temporal  mismatches  between  ther¬ 
mal  production  by  the  SOFC  micro-cogeneration  device  and  space 
heating  demands.  During  this  period,  the  make-up  burner  often 
cycles  on  in  the  early  morning  to  meet  space  heating  demands, 
while  the  heat  dumping  facility  is  activated  later  in  the  day.  Fully 
36.5%  of  the  SOFC  micro-cogeneration  device’s  thermal  output  is 
wasted  through  the  heat  dump  facility  while  another  9%  is  lost  due 
to  heat  transfer  from  the  tank’s  external  surfaces  to  the  containing 
room. 
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Fig.  13.  Monthly  integrated  electrical  balance  for  house  with  SOFC  micro¬ 
cogeneration  device. 


These  simulation  results  identify  a  number  of  areas  where 
the  performance  of  the  SOFC  micro-cogeneration-based  heating 
system  might  be  improved  to  increase  overall  system  efficiency. 
Further  analysis  with  the  calibrated  model  could  lead  to  improve¬ 
ments  in  the  design  and  operation  of  the  coupling  between  the 
SOFC  micro-cogeneration  device  and  the  heating  plant.  Alternate 
schemes  for  utilizing  the  SOFC  micro-cogeneration  device’s  ther¬ 
mal  output,  such  as  DFIW  heating  and/or  thermally  activated 
cooling,  could  be  explored.  The  potential  impact  of  design  improve¬ 
ments  to  sub-systems,  such  as  the  DC-AC  power  conditioner  or 
gas-to-water  heat  exchanger,  could  be  explored  by  adjusting  some 
of  the  calibration  coefficients  determined  in  Section  5.  Hypothetical 
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Fig.  14.  Monthly  integrated  thermal  energy  balance  on  tank  coupled  to  SOFC  micro¬ 
cogeneration  device. 


devices  could  be  analyzed  if  methods  existed  to  establish  appro¬ 
priate  coefficients  for  the  model’s  sub-systems.  This  would  enable 
analyses  focused  upon  optimizing  the  output  capacity  and  control 
of  SOFC  micro-cogeneration  devices  at  early  stages  of  design. 

7.  Conclusions 

Previous  research  has  produced  a  mathematical  model  for 
simulating  the  thermal  and  electrical  performance  of  SOFC 
micro-cogeneration  devices  and  which  is  suitable  for  use  in  whole¬ 
building  simulation  programs.  This  is  a  system-level  model  that 
considers  the  thermodynamic  performance  of  all  components 
that  consume  energy  and  produce  the  SOFC  micro-cogeneration 
device’s  thermal  and  electrical  output.  Each  of  the  model’s  nine 
control  volumes  represents  a  sub-system  that  produces  electrical 
power,  supplies  air,  captures  heat  from  the  hot  product  gases,  etc. 

Through  example,  this  paper  has  demonstrated  how  this 
model  can  be  calibrated  to  represent  the  performance  of  specific 
SOFC  micro-cogeneration  devices.  The  experimental  procedures 
that  were  employed  to  test  a  prototype  2.8kWAC  SOFC  micro¬ 
cogeneration  system  were  described  in  detail.  The  experimental 
configuration,  types  of  instrumentation  employed,  and  the  oper¬ 
ating  scenarios  examined  were  treated.  The  propagation  of 
measurement  uncertainty  into  the  derived  quantities  that  are 
necessary  for  model  calibration  was  demonstrated  by  focusing 
upon  the  SOFC  micro-cogeneration  system’s  gas-to-water  heat 
exchanger.  The  calibration  of  the  model’s  empirical  coefficients 
was  then  demonstrated.  The  results  of  this  calibration  approach 
were  presented  for  the  gas-to-water  heat  exchanger,  and  for  the 
fuel  cell  power  module’s  electrical  conversion  efficiency  and  air 
supply.  Although  not  presented  here,  the  same  methods  and  exper¬ 
imental  data  were  used  to  calibrate  the  other  aspects  of  the  model 
for  this  prototype  device,  namely  its  transient  response,  its  DC-AC 
power  conditioning  system,  its  dilution  air  system,  and  its  thermal 
losses.  The  reader  is  cautioned  that  the  calibrated  empirical  coef¬ 
ficients  presented  in  this  paper  are  only  valid  within  the  ranges  of 
independent  variables  examined  in  the  experiments. 

The  result  of  this  work  is  a  calibrated  model  that  can  accurately 
simulate  the  thermal  and  electrical  performance  of  this  2.8 1<WAC 
SOFC  micro-cogeneration  device  when  coupled  to  a  house’s  ther¬ 
mal  plant.  An  illustration  of  the  types  of  analyses  that  have  been 
enabled  with  this  work  was  given  by  focusing  upon  the  integra¬ 
tion  of  a  SOFC  micro-cogeneration  device  with  a  house’s  space 
heating  system.  This  demonstration  of  the  modelling  capabilities 
revealed  a  number  of  areas  where  the  performance  of  the  SOFC 
micro-cogeneration-based  heating  system  might  be  improved  to 
increase  overall  system  efficiency  and  described  how  the  model 
could  inform  the  design  and  development  process  for  SOFC  micro¬ 
cogeneration  devices. 

A  total  of  45  experiments  were  conducted  to  gather  the 
necessary  data  to  calibrate  the  model  for  the  2.8 1<WAC  SOFC  micro¬ 
cogeneration  device.  The  results  presented  here  have  demonstrated 
that  the  mathematical  form  of  the  individual  terms  of  the  model 
well  represents  the  functional  dependency  upon  the  independent 
variables.  However,  these  results  do  not  speak  to  the  validity  of  the 
calibrated  model  to  accurately  predict  performance  under  operat¬ 
ing  regimes  other  than  those  examined  during  the  45  experiments. 
In  fact  16  additional  experiments  were  conducted  (disjunct  from 
the  45  experiments  used  to  calibrate  the  model)  to  empirically  val¬ 
idate  the  model  and  its  calibration,  the  results  of  which  will  be 
reported  in  a  future  paper. 

Future  research  will  apply  the  calibrated  model  to  examine  the 
thermal  and  electrical  performance  of  the  2.8 1<WAC  SOFC  micro¬ 
cogeneration  device  in  various  residential  buildings,  in  various 
climates,  and  under  various  operating  scenarios  to  contrast  the  per- 
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formance  of  the  technology  to  conventional  heating  and  electrical 
supply  systems.  Using  the  techniques  described  here,  the  model 
will  also  be  calibrated  to  represent  the  performance  of  other  proto¬ 
type  SOFC  micro-cogeneration  devices  and  configured  to  represent 
hypothetical  devices  to  guide  research  and  development  efforts. 
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